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The geometry of optical lattices can be engineered allowing the study of atomic transport along paths ar¬ 
ranged in patterns that are otherwise difficult to probe in the solid state. A question feasible to atomic systems is 
related to the speed of propagation of matter-waves as a function of the lattice geometry. To address this issue, 
we have investigated theoretically the quantum transport of non-interacting and weakly-interacting ultracold 
fermionic atoms in several 2D optical lattice geometries. We find that the triangular lattice has a higher propa¬ 
gation velocity compared to the square lattice, and the cross-linked square lattice has an even faster propagation 
velocity. The increase results from the mixing of the momentum states which leads to different group velocities 
in quantum systems. Standard band theory provides an explanation and allows for a systematic way to search 
and design systems with controllable matter-wave propagation. Moreover, the presence of a flat band such as 
in a two-leg ladder geometry leads to a dynamical density discontinuity due to its localized atoms. Possible 
realizations of those dynamical phenomena are discussed. 
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I. INTRODUCTION 

Variational calculus leads to the conclusion that the shortest 
distance a free classical particle follows between two points is 
a geodesic HI. Is this also true for quantum particles? As 
expected the answer is not trivial. For instance, in the stud¬ 
ies of quantum dynamics of magnetic domains, it has been 
found that interesting patterns emerge 0, and an analogue 
of a quantum-mechanical “forest fire propagation” has been 
realized ||3l. Moreover, when the classical random-walk prob¬ 
lem is promoted to a quantum one, the dynamics is no longer 
diffusive and the object spreads faster (see iH and references 
therein). Therefore, quantum effects, in particular the wave 
nature, can alter the dynamics in a fundamental way. The de¬ 
tection of such unusual behavior can be, however, difficult in 
the solid state. To measure matter waves representing elec¬ 
trons in metals, attosecond spectroscopy needs to be imple¬ 
mented and the observation length scale is limited by scat¬ 
tering effects Q . Recent demonstrations of controlling group 
velocity of cold-atom matter waves ||6l and photons 0 further 
illustrate the feasibility of tuning group velocity in quantum 
systems. 

Ultracold atoms in engineered optical lattices allow direct 
measurement of matter-wave dynamics. Many interesting 
2D lattice geometries, including square, triangle, honeycomb, 
kagome, have been fabricated lIS UT^ . These systems have 
been shown to be versatile quantum simulators of complex 
many-body systems lfT3UT6ll because of the wide selections of 
atomic species and tunable parameters such as interactions or 
trapping potentials. Recent advances in introducing time de¬ 
pendence to the parameters and configurations of cold-atom 
systems open up opportunities for studying nonequilibrium 
physics, particularly transport phenomena inHii]. Impor¬ 
tantly, while the Fermi velocity, vp, of electrons in typical 
metals such as copper is on the order of 10® m/s, in cold-atom 
systems vp is on the order of 10“® m/s ll2^ . Such a slow 
motion of cold atoms then allows detailed analyses of their 


dynamics. 

As illustrated by the matter-wave example, transport of 
fermionic atoms is closely related to electronic transport in 
nanoscale and mesoscopic systems ||20l|2ll|23llMl- For in¬ 
stance, when ultracold fermions are driven out of equilibrium 
in a one-dimensional optical lattice, a quasi-steady state cur¬ 
rent with a constant magnitude for a period of time emerges, 
which is the precursor of the steady-state current (in the ther¬ 
modynamic limit) found in biased solid-state systems ll^l^ . 
Quasi-steady-state currents have also been found to survive in 
two-dimensional systems Eli. Interestingly, this quasi-steady 
state persists also in the absence of particle interactions, a fact 
not easily verifiable in the solid state, while the presence of 
strong interactions can change the transport from ballistic to 
diffusive E^ . 

Simulating band-theory of electrons using cold-atoms has 
made it possible to unambiguously demonstrate paradigmatic 
phenomena. For example, Bloch oscillations are known in 
band theory but defied a clear demonstration in solid state sys¬ 
tems until their demonstrations in cold-atom systems (see ll29l 
and references therein). It is thus reasonable to expect that 
cold-atom systems can provide more insights into geometric 
effects on matter-wave propagation. Since a mass current of 
cold atoms corresponds to a traveling matter wave, it is inter¬ 
esting to clarify how fast the wavefront propagates, and how 
its speed can be controlled, thus providing an answer to the 
question posited at the beginning. 

To explore the dependence of group velocity on the un¬ 
derlying geometry, we have employed noninteracting and 
weakly-interacting fermions in various optical lattice geome¬ 
tries and found that quantum matter waves can be accelerated 
altering the underlying lattice geometry, which is not expected 
in classical particles. Using band theory in the thermodynamic 
limit, we found that interference of matter waves plays a key 
role in the atom dynamics. This also provides guidance for 
finding lattice geometries with faster matter-wave propaga¬ 
tion. Including a weak repulsive interaction at the mean-field 
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level does not change the conclusions qualitatively. 

We also note that certain lattice geometries can support flat 
bands, which refer to a special class of dispersionless bands 
(see Ref ||30H^ and references therein). For example, the 
kagome lattice can support a flat band and can be realized in 
optical lattices iia. The particles residing on a flat band do 
not possess kinetic energy and as a consequence, they do not 
participate directly to transport and lead to interesting phe¬ 
nomena. For instance, a dynamically generated flat-band in¬ 
sulator sustaining a density discontinuity has been predicted in 
optical kagome lattices Ezl. We find transport features simi¬ 
lar to this flat band effect in a two-leg ladder, and discuss its 
possible experimental realization. Such examples suggest that 
not only the mobile properties but also interesting insulating 
ones could be explored using cold atoms in engineered lattice 
potentials. 

In addition, the zig-zag lattice provides an opportunity to 
explore rich physics related to frustrations in optical lattice 
systems 133] [Mj. One important feature of the zig-zag lattice 
is that the roles of the nearest neighbors and the next nearest 
neighbors can be switched by tuning the coupling or lattice 
parameters. We explore the same idea in transport by consid¬ 
ering noninteracting fermions in a zig-zag lattice with tunable 
tunneling coefficients and show that the speed of matter-wave 
propagation can be controlled. A maximal velocity emerges 
as the relative tunneling coefficients are tuned continuously. 

This paper is organized as follows. In Sec. |n| the lat¬ 
tice geometries and models for studying matter-wave prop¬ 
agation are introduced, and band theory is briefly reviewed. 
We present the geometry-dependent matter-wave velocity and 
flat-band induced density discontinuity, along with their ex¬ 
planations from band theory in the thermodynamic limit, in 
Sec. Ill Sec. IV discusses possible experimental setups for 
verifying the results and applications in controlling matter- 
wave propagation. Then Sec. |V| concludes this work. 


II. THEORETICAL BACKGROUND 

To highlight interesting geometrical effects on quantum 
transport, we consider three types of 2D lattices illustrated in 
Figure They are the square lattice, triangular lattice, and 
square lattice with next-nearest neighbor hopping (the SNNN 
lattice). The SNNN lattice can be viewed as a Bravais lattice 
with two lattice sites per unit cell, which we may call A and 
B. In the SNNN lattice both the nearest neighbor (A-A and 
B-B) and next-nearest neighbor (A-B) hoppings are allowed. 
One may consider the SNNN lattice as the 2D version of the 
body-centered cubic lattice in 3D llTSll because the SNNN lat¬ 
tice can be constructed by two intercalating square sub-lattices 
(labeled as A and B). The SNNN lattice is also geometrically 
similar to the checkerboard lattice demonstrated in Ref. im. 
In our discussion the relative hopping coefficients of the A-A, 
A-B, and B-B links are assumed to be tunable. A possible 
experimental realization may use bilayered square lattices as 
shown in Fig.[TJd), although the need for a phase shift between 
the two layers could be a challenge. By viewing perpendicu¬ 
larly into the bilayer lattice, the two square sublattices (A and 


B) are intercalated as depicted by the SNNN lattice explained 
in Figure [T](e) and (f). The lattice constant a serving as the 
length unit is chosen to be the same for the square and trian¬ 
gular lattice and the A-B link in the SNNN lattice. 

For noninteracting single-component fermions in a moder¬ 
ate lattice potential, the system may be modeled by a tight- 
binding Hamiltonian of the form ll36l 

H = — iijc\cj, (1) 

(U> 

where {ij) denotes a pair of sites connected by a link, Ci (cl) 
annihilates (creates) a fermion at site i, and tij is the hop¬ 
ping coefficient. For a uniform lattice with tij = t, the unit 
of time is defined as to = h/i. We set h = 1 and assume 
that there are (Ly) lattice sites along the horizontal (verti¬ 
cal) direction. Recently developed box potentials ifJTl make it 
practical to study homogeneous properties in cold-atom sys¬ 
tems. Moreover, a weak background harmonic potential does 
not change transport properties qualitatively ||25]| . Therefore, 
we will focus on the intrinsic transport phenomena in a homo¬ 
geneous system with open boundary conditions and employ 
the microcanonical picture of transport which is ideally suited 
for closed finite systems ll^lMl . 

The initial condition is similar to that in Ref ll26l . where 
the system is separated into two regions (left and right) by 
a laser sheet, and fermions are only loaded to the left, as il¬ 
lustrated in Fig. We first consider single-species fermionic 
atoms, which are noninteracting fermions. In Refs, 
it has been shown that the dynamics is quite insensitive to 
the initial filling as long as the corresponding thermodynam¬ 
ics limit exists, so we consider an initial band insulator on the 
left half lattice. Adjusting the initial filling should not lead 
to qualitative differences, except in the flat-band effect dis¬ 
cussed later. At time t = 0 the optical barrier is lifted and 
the fermions start to propagate to the right. The full quan¬ 
tum dynamics of noninteracting fermions can be monitored 
by the single particle correlation matrix defined by its ele¬ 
ments Cij{t) = {cl{t)cj{t)) ll25l and from this we compute 
the density on site i, ni(t) = Cii{t). 

The equations of motion for (t) can be obtained via 

.9(C^ = {[clcj,H]) = ^(fi,A(4_ACj)-ij,A(clcj+A)). 

A 

( 2 ) 

Here A denotes the vector to the other side of a link. The 
equation of motion was evaluated numerically using the 
fourth-order Runge-Kutta method The initial correla¬ 

tion matrix was set up with one particle located at each site on 
the left half of the lattice. We have calibrated the numerical 
procedure using known exact solutions. 

In the thermodynamic limit of infinitely large lattices, the 
results should agree with band theory ll35l . By implementing 
the lattice version of the Fourier transform, the Hamiltonian 
becomes 


H = '^ekclck. (3) 

k 
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Figure 1. (Color online) Illustration of the lattice geometries: (a) 
Two-leg square lattice. The coordinates x, y are also shown, (b) 
Two-leg triangular lattice, (c) Two-leg square lattice with next- 
nearest neighbor hopping (SNNN) lattice with the sub-lattices A and 
B labeled. The hoppings between the A-A, A-B, and B-B links are 
all allowed in the SNNN lattice. Initially the left half of the lattice is 
filled (solid circles) and the right half is empty (empty circles). The 
density imbalance drives a matter-wave propagating to the right-hand 
side, indicated by the velocity vector, (d) Possible experimental real¬ 
ization of the SNNN lattice with tunable tunnelings using bilayered 
optical lattices. One layer has a square lattice of the A-sites and the 
other layer has a square lattice of the B-sites. When viewed from the 
top of the bilayer lattice, inter- and intra- layer tunneling constmcts 
the A-A, A-B, and B-B links as shown in (e). (f) By rotating the 
geometry of (e) by tt/T and using two laser sheets (denoted by the 
thick red lines) to confine a two-leg region, the SNNN lattice of (c) 
can be realized. 


Here Cfc (c^) is the annihilation (creation) operator in momen¬ 
tum space. The dispersion can be exactly solved for non¬ 
interacting fermions. Importantly, the semiclassical group ve¬ 
locity is given by llTSll 


to how each site is connected to its neighbors, with suitably 
chosen lattice geometry the propagation speed can be altered. 


Figure shows the evolution of the y-direction averaged 
density, n{ix) = '^i^,iy)lLy, for the square, triangu¬ 

lar, and SNNN lattices with {Lx,Ly) = (200, 2) ((a)-(c)) and 
for the SNNN lattice with {Lg;,Ly) = (200,3). At f = 0 
the left half is filled and the right half is empty. The bal¬ 
listic transport of fermions clearly shows a light-cone struc¬ 
ture Eomii with the outermost wavefront showing a con¬ 
stant maximal propogating speed. For different geometries 
the maximal speeds are different, and there are equivalent 
ways for measuring them. For instance, the expansion rate 
of the perturbed region of the cloud given by Aw/At will 
be twice of the maximal propagating speed m and can be 
extracted following the illustration of Fig. (e). Since the 
ballistic expansion of noninteracting fermions leads to con¬ 
stant maximal propagating speed as shown in Fig. 0 (a)-(d), 
measuring the time tx it takes for the outermost wavefront 
to reach a point ix can also be used to evaluate the maximal 
speed Vm = {ix — iL/ 2 )/tx- In our simulations we impose 
open boundaries and use the latter method to determine the 
maximal propagating speed. A density-averaged velocity may 
be defined following Ref. ED, which works as an indicator 
of the average velocity of condensed bosons. Since fermions 
span the momentum space up to the Fermi momentum, the 
density-averaged velocity is less informative in revealing the 
spreading velocity and here we focus on the maximal group 
velocity of fermionic matter waves. 


It can be observed that the wavefront speed increases as 
more diagonal links are added to the lattice. The full quantum 
dynamics thus claims that altering the lattice geometry may 
boost matter-wave propagation, which can be very difficult to 
explain using classical physics. Even more puzzling is that 
in the case of the SNNN lattice with Lx = 200 and Ly = 
2, about half of the particles stay on the left half, and this 
causes a density discontinuity at Lx 12, which is visible on 
Figure l^c). In contrast, for the SNNN lattice with Lx = 200 
and Lj, = 3 no density discontinuity can be observed. This 
indicates that Lj, = 2 is special for the SNNN lattice. 


v/c = Vfcefe. (4) 

This semiclassical prediction in the thermodynamic limit will 
be compared to the fully quantum mechanical results of 
matter-wave propagation. 

III. RESULT AND DISCUSSION 

We first consider uniform lattices where all hopping coef¬ 
ficients are equal: tij = i. The time it takes for a particle 
to tunnel through an isolated double-well potential is 7rfo/2 
ll40l . and in a ID multi-site lattice it asymptotically takes fo/2 
to tunnel through one link ll25l l26l because of the modifica¬ 
tion of the dispersion. Since quantum matter wave is sensitive 


By analyzing ri as a function of t as shown in Fig. one 
can estimate the time it takes for the wavefront to first arrive 
at a point ix- For convenience, we take ix = Lx and measure 
t* as the time when ri reaches a threshold and extract the ve¬ 
locity V = {Lx/2)/t* for the three types of lattices shown in 
Fig .[^as a function of Ly with Lx = 200 and present the data 
in Figure[^(a). The details of ri as a function of t for selected 
cases are shown in Figure|^(bl) and (b2). The threshold is set 
to be 0.05 and the results do not change qualitatively with this 
choice. We remark that by measuring the spreading rate of 
the density profile shown in Fig. and using the same thresh¬ 
old to determine the boundary of the density spreading, the 
same maximal propogating speed would be determined. The 
agreement is guaranteed by the light-cone spreading of nonin¬ 
teracting fermions shown in Fig. 
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Figure 3. (a) Matter-wave propagation velocities for the square (cir¬ 
cular symbols), triangular (square symbols), and SNNN (diamond 
symbols) lattices for = 200 and selected Ly. The horizontal lines 
show the maximal velocities calculated from band theory. The evo¬ 
lutions of the j/-direction averaged density n for two selected cases 
indicated as b. 1 (for the SNNN lattice) and b.2 (for the square lattice) 
are shown in panels b. 1 and b.2. The wavefront arrival time at 
(denoted by t*) is determined by the time when h rises above a se¬ 
lected threshold. In this work we choose the threshold to be 0.05 as 
shown by the dashed lines in panels b.l and b.2. 
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A. Band theory 

The geometry-dependent propagation velocity can be ac¬ 
counted for by the standard band theory for fermions llTSl . 
For an infinitely large 2D uniform square lattice following the 
illustration of Fig. (a), we introduce lattice Fourier trans¬ 
form with lattice momentum k = {kx,ky) and diagonalize 
the Flamiltonian ll29l . Then the energy dispersion is 

es(^) = —2f[cos(fca;) + cos(A:y)]. (5) 

The semiclassical velocity in the a;-direction is 

uf = 2tsiii{kx)- (6) 

Thus the maximal velocity in the x-direction is 2t, which is 
consistent with the numerical simulations for the square lat¬ 
tice. Since kx and ky are decoupled, the result is independent 
of ky so the maximal velocity already reaches the band-theory 
prediction when Ly = 2, as shown in Fig.[^(a). 

For the triangular lattice, the energy band is 


Figure 2. (Color online) Evolution of the particle density. Here the 
brighter regions correspond to higher density, (a) Square lattice, (b) 
Triangular lattice, (c) Two-leg SNNN lattice. Here Lx = 200 and 
Ly = 2. In (c) there is a density discontinuity at Lxl2 caused by the 
fiat band, (d) The SNNN lattice with Lx = 200 and Ly = 3 shows 
no density discontinuity, (e) Using (a) as an example, the spreading 
of density, Aui/ At, gives twice the maximal group velocity. 


trik) = —2t 

The corresponding velocity is 


n \ , r, { kx x / '\/3ky . 

cos[kx) + 2cos(y)cos( ^ 


uj = 2t 


• / 7 ^ I ' / \ / “'^^ky , 

sm[kx) + sin(y) cos(^-^) 


(7) 


( 8 ) 
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Unlike the square lattice, however, an important feature for 
this case is that kx and ky are mixed in the expressions. The 
maximal velocity is 2f(l + \/2j2), which is larger than that 
in the square lattice. This ky dependence also suggests that 
the Ly dependence will be more significant. Indeed, Fig.|^(a) 
shows a stronger Ly dependence of the wavefront propagation 
velocity for the triangular lattice. 

For the SNNN lattice, we hrst consider the inhnite 2D lat¬ 
tice and defer the discussion of the special case of Ly = 2 to 
later. Since there are two intercalated sublattices (A and B), 
the energy band will split into two bands while the size of the 
first Brillouin zone is only half when compared to the previous 
two cases. The energy dispersion is 

esNNN{k) = -‘ltcos{kx) ± tg{k). (9) 

Here g{k) = a/1 + 4 cos{kx) cos{ky) + 4:COs‘^{kx)- The 

group velocity is 

ySNNN _ cos{ky) =F 

t[2sin(fca;) cos{ky) + icos^kx) sm{kx)]/g{k). 

( 10 ) 

A plot of the two bands at ky = n/a is shown in Figure 
The two bands touch at four points in the first Brillouin zone 
so the system remains conducting for any filling less than the 
band insulator. The maximal velocity for the SNNN lattice is 
6t. 

As shown in Fig. the wavefront velocity of the SNNN 
lattice approaches the maximal velocity as Ly increases, but 
the convergence is slower when compared to the previous two 
cases. Fig. |^d) shows many trailing wakes behind the first 
wavefront in the SNNN lattice case and they are the byprod¬ 
ucts of the more complex dispersion 0. As a consequence, 
the averaged density n of the SNNN lattice rises less abruptly 
when compared to the other cases as shown in Fig. [^(b.l) 
and (b.2). By estimating the velocity from the time when 
n rises above the selected threshold 0.05, an underestima¬ 
tion is in place. This can be improved by choosing a lower 
threshold of n, although it would require more precise exper¬ 
iments to locate the lower-density threshold. As illustrated 
in Fig. 1^ the same issue may persist if one measures the 
density-averaged velocity ED. The low-density wavefront 
corresponds to larger distances from the initial boundary, so 
the product of the density and the distance around the wave- 
front leads to higher uncertainties in the averaged velocity. 

We have two remarks. Our initial condition corresponds to 
a band insulator on the left half lattice. It has been shown 
that reducing the initial filling only reduces the magnitude of 
the current and the dynamics remains similar E5l WH . One 
may use an initial configuration with finite total momentum, 
which may be realizable by filling up fermions in selected 
momentum states or imposing a momentum transfer. Since 
the spectral weight of the states are determined by the initial 
condition, the propagation will be biased towards the selected 
initial momentum. The maximal propagating speed, however, 
is determined by the lattice geometry and will not be changed 
by different initial states. Here we focus on fermionic trans¬ 
port because it has been shown that noninteracting bosons in 



ka ka 

Figure 4. Energy bands of (a) 2D SNNN lattice with ky = nja and 
(b) two-leg SNNN lattice where a flat band emerges from the over¬ 
lapping of energy bands when the tunneling coefficients are uniform. 


the ground state will form a macroscopically coherent conden¬ 
sate and the transport is highly oscillatory ll25l . In addition, 
at low temperature bosons stay in the bottom of energy bands, 
so they are insensitive to flat-bands which are at higher en¬ 
ergies. For the purpose of probing the maximal propagating 
speed and flat-band physics, fermions are then a better choice. 

The honeycomb lattice, when the zig-zag side is placed 
along the horizontal direction, can be viewed as a deformed 
square lattice with alternating missing vertical links. In this 
case, the wavefront propagation velocity from the band the¬ 
ory and simulations are both identical to that of the square 
lattice when geometrical factors are properly considered. It is 
known that the low-energy dispersion in the honeycomb lat¬ 
tice around the Dirac points resembles relativistic fermions 
E^ . For global transport studied here, the low-energy effects 
are not observable. 


B. Flat band effects 


We now turn to the reason behind the density jump of the 
two-leg SNNN lattice observed in Fig. When Ly = 2, the 
Hamiltonian can also be diagonalized and the energy bands 
are 


i{kx) = —2tcos{kx) ± t\l -F 2cos(fca;)|. (11) 


The two bands are shown in Fig. j^b). One important fea¬ 
ture is that there is a flat band, i.e., a collection of disper¬ 
sionless states. The particles on the flat band have localized 
spatial patterns so they do not participate in transport directly. 
For the two-leg SNNN lattice, the flat-band state corresponds 


to the wavefunction ^ ^ localized on a vertical link. 

The tunneling from the two sites on a verticle link to any site 
to the left or to the right leads to a total destructive interfer¬ 
ence, which explains the localization. This also contrasts why 
Ny = 2 is special for the SNNN lattice because no such total 
desctructive interference can be found in other values of Ny. 
In the presence of a chemical potential difference caused by 
the initial density difference, mobile atoms on the curved part 
of the dispersions are driven to the initial vacuum region. 

The atoms on the flat-band, in contrast, remain in the ini¬ 
tially filled region because their wave functions tend to lo- 
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calize. Since each momentum state can accommodate one 
fermion, in Fig. |^b) the flat-band contains exactly states. 
Thus, the flat-band states account for half of the total states 
(i.e.. Lx out of the total 2Lx states). Therefore, the static 
property of the localized flat-band states is the reason for the 
density jump. By inspecting the particle density profiles at dif¬ 
ferent time slots, we have verified that the density jump does 
not decrease significantly as the wavefront propagates to the 
right. We also remark that, because fermions pile up from the 
lowest energy to the Fermi energy, the lower part of the dis¬ 
persive band will be populated first. There are {2/3)Lx states 
on the dispersive band below the energy of the flat-band as 
shown in Fig.Qb). Thus, to observe the density discontinuity 
the initial Ailing on the left half lattice should exceed 1 /3 for 
the flat-band to be filled. This analysis again shows the power 
of band theory in explaining transport properties. 

Such a flat-band induced dynamical density discontinu¬ 
ity may be generic in cold-atom systems. For instance, the 
kagome lattice has one flat band in its lowest three bands and 
quantum dynamics show that a density jump of magnitude 1 /3 
emerges as particles flow from an initially filled region into an 
initially empty region There is, however, a subtle differ¬ 
ence between the flat band of the kagome lattice and that of 
the two-leg SNNN lattice. The kagome flat-band is on top of 
the other two dispersive bands while that in the two-leg SNNN 
intersects the other dispersive band. Nevertheless, their roles 
in transport are identical. While demonstrating such a dynam¬ 
ical density discontinuity is straightforward in cold-atoms sys¬ 
tems, in conventional condensed matter depleting mobile elec¬ 
trons completely can be very challenging and observing this 
phenomenon can be a daunting task. 

A few remarks on the two-leg SNNN lattice are worth men¬ 
tioning. First, the existence of the flat band depends crucially 
on the condition of uniform tunneling coefficients. If the di¬ 
agonal (A-A and B-B) links have a different tunneling coeffi¬ 
cient F, the bands start to curve and there is no longer a flat 
band. We have tested the case for F < i and found that when 
F/t is close to 1, the flat-band is bent only slightly and a sharp 
density difference is still observable, but as F/f decreases, the 
sharp density difference decays away. The resulting disper¬ 
sion of the two-leg SNNN lattice 

e(k) =—2icos{kx) ± \i + 2F cos(fca:)| (12) 

yields two maximal group velocities for the two bands. As the 
flat band distorts, the atoms propagate in a light-cone structure 
bound by v+ = 2{t — F), while atoms in the dispersive band 
reach a maximum velocity of z;_ = 2{t + F). Figure [^sum¬ 
marizes the behavior as F/i is decreased. We also caution that 
the density-averaged velocity ED may not fully account for 
the underlying transport when multiple bands are present. For 
instance, there are two maximal velocities shown in Fig. [^ 
and an averaged velocity will average out this feature. 

Secondly, in the presence of magnetic flux penetrating the 
ladder, the dispersion exhibits additional interesting features 
and the system is known as the Creutz ladder ll42l l43l and 
transport in the square lattice with effective magnetic flux has 
been recently studied in cold-atom systems ll44l . Here, we 
focus only on geometrical effects on transport. 



site 


Figure 5. Density evolution of two-leg SNNN lattices with F/i = 
0.9 (a), 0.8 (b), and 0.7 (c). Here Lx = 200. The case with F/t = 
1 has been shown in Fig. |^(c). (d) Using (b) as an example, the 
maximal rates of density spreading, Aw±/At±, bound the light- 
cone structures of the two bands. 


C. Weak interactions 

Transport behavior can be altered by introducing strong in¬ 
teractions as demonstrated in Ref ESlI . In the presence of in¬ 
teraction, however, band theory has limited applicability and 
full numerical simulations are usually required. A lack of 
efficient numerical methods for studying 2D interacting sys¬ 
tems makes a full discussion of interaction effects out of the 
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scope of this work. Here we only demonstrate that in the pres¬ 
ence of weak interactions which may be treated in the mean- 
field approximation, the geometrical effects are still observ¬ 
able. We consider two-component fermions labeled by spins 
cr =t, i, which may be different species of atoms or the same 
atoms in two different internal states. For contact interactions 
common in cold atoms, the system may be described by the 
Hubbard model l36l + where 

nia = cl^Cia and U is the onsite repulsion coupling constant. 

The equations of motion for the Hubbard model can¬ 
not be solved analytically due to the presence of multi¬ 
particle correlations. Here, we follow Ref. Il26l and imple¬ 
ment the standard Hartree-Fock approximation by decom¬ 
posing (c|^^Ci^c|'^C 4 ) as {c\^Ci-\){c\^Ci]). Then the equa¬ 
tions of motion become idt{clcj) = X)A(^i.A(c-Cj+A) + 
U,A{4-A'^j)) + Ui{4cj){4c^) - {4cj){c]cj)), which can 
be solved numerically. We impose the condition = 

{4iCii) at any time. Our numerical results for 0 < 17 < t 
show no qualitative difference from the corresponding nonin¬ 
teracting fermion cases. On the other hand, the interactions 
could influence the flat band and there will be leaks from the 
flat-band states. However, the leaking rate is extremely small 
in the weakly interacting regime U/i < 1. Within the time 
scale of our simulations, the density discontinuity remains vis¬ 
ible. Therefore, the conclusions from the noninteracting case 
should hold when the background interactions are weak. 


IV. EXPERIMENTAL IMPLICATIONS 
A. Possible realization of the SNNN lattice 

Here we discuss possible experimental realizations and ap¬ 
plications of this work. The optical square and triangular 
lattices have been demonstrated EHia and geometrical ef¬ 
fects on transport should be observable on those lattices. The 
SNNN lattice is similar to the checkerboard lattice of Ref. im 
whose diagonal tunneling (A-A or B-B) is expected to be 
weaker than the nearest-neighbor tunneling (A-B) due to the 
in-plane geometry. The SNNN lattice with a tunable ratio be¬ 
tween those two tunneling coefficients may be created by bi¬ 
layer square lattices illustrated in Figure[2(e), which is similar 
to those discussed in Ref. 1451 . but the need of an additional 
phase shift between the two layers will be a challenge. By 
adjusting either the lattice depth on each sheet or the distance 
between the two sheets, the relative strength of the two tun¬ 
neling coefficients can be tuned. 

While genuine bilayer optical lattices have not be available, 
multi-layer optical lattices have been experimentally realized 
for the honeycomb lattice as an analogue of graphene ll46l . 
To realize two-leg ladders of the SNNN lattice for supporting 
the flat-band, a tight background harmonic trap or laser sheets 
Ell for confining the atoms in a narrow transverse region may 
also be needed. As the geometry and tunneling coefficients 
change, the energy dispersion as well as the group velocity of 
the fermions change accordingly. 



Figure 6. (Color online) Matter-wave propagation velocity of a zig¬ 
zag lattice (illustrated in the inset) with tunable tunneling coeffi¬ 
cients. The zig-zag lattice has two sets of tunneling coefficients: i 
for the horizontal tunneling and t' for the diagonal tunneling. When 
t' is fixed, the velocity increases as the ratio f/T (circle symbols) 
increases from 0 to 1. When t is fixed, the velocity increases 2 &t' jt 
(triangular symbols) increases from 0 to 1. A maximal group veloc¬ 
ity emerges when t' jt = 1. 


B. Possible application 

We end our discussion of geometrical effects on matter- 
wave propagation by presenting a possible application in the 
zig-zag lattice optical lattice ||33]|34l, which is basically a two- 
leg ladder of the triangular lattice as shown in the inset of Fig¬ 
ure]^ The hopping coefficient along the horizontal direction 
is t and that along the diagonal is F, and the ratio between 
them are assumed to be continuously tunable. The initial con¬ 
dition is identical to the previous cases with the left half of the 
lattice being a band insulator and the right half being empty, 
as illustrated in Fig. 

When F = 0 and t gives the time scale the two legs 
decouple and it takes for the matter wave to reach the 

far right boundary. In the opposite limit when 7=0 and 
t' gives the same time scale it takes Lxto/2 because the 
particle has to traverse twice the distance to reach the far right 
boundary. The velocities are thus 2a /to and a /to for the two 
limits. Figure|^shows the velocity for selected values of t/t' 
and i'/i. We chose the quantity in the denominator to give the 
same tunneling time scale to to guarantee a fair comparison. 
Importantly, a maximal velocity appears when t = F when 
each site has the most available number of neighbors to tunnel 
into. 

By tuning the ratio of the tunneling coefficients, a zig-zag 
optical lattice may serve as an atomic “bicycle gearing” for 
changing the speed of atoms passing through it. One pos¬ 
sible scenario is to control the nearest-neighbor tunneling i 
while keep the next-nearest-neighbor tunneling F fixed. When 
7=0, the velocity is slowed down to a/to- To accelerate the 
matter wave, one can tune to 7 = 7' and the velocity will be 
increased by 270%, according to Fig. The zig-zag lattice 
tunes matter-wave propogating speed by the relative hopping 
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strength of the two types of links (horizontal and diagonal) 
while the total current through the lattice is roughly main¬ 
tained in a finite range because the total bandwidth does not 
change drastically. One may consider an alternative way of 
simply reducing all the hopping coefficients of a lattice to de¬ 
crease the overall propagating speed, but this will reduce the 
total current due to the reduced bandwidth as well. An ex¬ 
perimental set up proposed for the zig zag lattice by superpo¬ 
sitioning a triangular lattice and a strong optical superlattice 
miiTiiisi could make such a gearing device feasible. 

Similar applications may be relevant to quantum informa¬ 
tion transfer by controlling the group velocity of information 
carriers played by massive quantum particles. There have 
been experiments demonstrating information exchanges be¬ 
tween photons and cold-atom matter-waves ia, where the 
propagtion speed of the matter wave affects the overall trans¬ 
fer of the information. Our investigations into geometrical 
control of matter-wave propagation could provide addtional 
manipulations of quantum systems. 

V. CONCLUSIONS 

We have explored geometrical effects on quantum trans¬ 
port of cold atoms. With different energy dispersions from 
various lattice geometries, significantly increases of matter- 


wave propagation should be observable. With the assistance 
of band theory, searching and implementing lattice geometries 
for controlling transport can be performed efficiently. More¬ 
over, insulating phases sustaining a density difference is a 
generic feature for lattice geometries supporting a flat band, 
and the different roles of mobile and localized atoms in quan¬ 
tum transport have been elucidated. Similar studies can be 
performed in the recently realized optical Lieb lattice ll50l and 
the proposed optical saw-tooth lattice ll34ll . which also support 
flat-bands. 

Optimal designs of lattice geometries for controlling and 
regulating quantum transport of cold atoms could contribute 
useful elements to the thriving field of atomtronics ifSTl l52ll 
for simulating or complementing electronics. The feasibility 
of tuning matter-wave propagation using geometrical effects 
may also find applications in quantum quench dynamics uni 
|53l, where the relaxation rate after a quench may be controlled 
by geometrical effects. 
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